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J!^ : Abstract 

*> • It is shown that a single, strongly scattering impurity produces a bound or a 

CN 

^•^ ■ virtual bound quasiparticle state inside the gap in a d-wave superconductor. 

!>' 

f^ ■ The explicit form of the bound state wave function is found to decay exponen- 

^: 

^*^ ■ tially with angle-dependent range. These states provide a natural explanation 






O 
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of the second Cu NMR rate arising from the sites close to Zn impurities in 



'^ ■ the cuprates. Finally, for finite concentration of impurities in a d-wave super- 

a: 

O ■ conductor, we reexamine the growth of these states into an impurity band, 



and discuss the Mott criterion for this band. 



Effects of impurities on tlie properties of superconductors liave been investigated in great 
detail for low-temperature [|I| , heavy- fermion , and liigh-temperature superconductors [Q . 
Tlie main reason for the interest in the effects of impurities on the superconducting state is 
the fact that the superconducting properties are quahtatively modified by impurity atoms, 
depending, for example, whether they are magnetic or non-magnetic. In principle, this 
observation can be useful as a method of identifying the nature of the pairing state in 
superconductors. For example, any magnetic impurity will be a strong pair-breaker for (s- 
wave, (i-wave, etc.) spin-singlet superconductors, in accord with the generalized Anderson 
theorem. On the other hand, even scalar (non-magnetic) impurities are pair-breakers for 
"higher-orbital-momentum" states, such as a (i-wave pairing state. 

The two main approaches in understanding the effects of impurities in conventional {s- 
wave) superconductors rely either on the strong- or on the weak-scattering limit, (a) The 
Abrikosov-Gor'kov (AG hereafter) formalism [Q treats impurities in the Born approximation. 
Any impurity problem is characterized by two physical parameters: the phase shift 6q due 
to impurity scattering (which we assume to be s-wave) and the concentration of impurities 
nimp. In the AG approach, the only parameter entering the formalism is the scattering rate, 
T^^ = (2nijnp/7rA^o) sin^ 60, proportional to the product of the concentration and sin^ 60. Here 
A'"o is the normal-state density of states at the Fermi energy, and So = NqV is the s-wave 
phase shift for a weak impurity potential V. Therefore, in the limit of dilute concentration 
of strong magnetic impurities, the AG approach will yield a small average scattering rate 
1^. On the other hand, the local properties of the superconductor near an impurity site, 
such as the local density of states and the gap amplitude, will be modified dramatically. In 
this limit (b) the Yu-Shiba approach |^,^ should be used, which treats magnetic impurities 
in the unitary-scattering limit with the s-wave phase shift 60 ~ 7r/2. It was shown by Yu 
and Shiba that, in the unitary limit, a localized magnetic impurity, interacting with the 
spin density of conduction electrons at the impurity site, produces a true bound state inside 
the energy gap, |ti;| < Aq, where the density of states vanishes. Note that, in general, the 
overlap with the particle-hole continuum only allows virtual states to be formed with finite 



lifetime. The relation between this approach and the AG treatment was established in [|^, 
where it was shown that in the Born limit one recovers the AG results, and the bound state 
is indistinguishably close to the band edge. 

Our work is partially motivated by the fact that non-magnetic impurities are known to 
be strong pair-breakers in a nontrivial superconductor. The Zn substitutions in cuprates are 
one example of this. Although Zn ions are nominally non-magnetic, Tc is strongly suppressed 
by Zn substitution of Cu in the planes P,^]. Therefore, it is reasonable to assume that Zn 
ions are non-magnetic unitary scatterers. 

The purpose of this Letter is to address the question of virtual impurity-bound states 
in a (i-wave superconductor and, within this framework, to explore possible implications of 
the assumption that the pairing in cuprates is in the dr^2_y2 channel. We shall generalize the 
original Yu-Shiba [^|^ approach to the case of arbitrary-strength of non-magnetic impurities 
in a da;2_y2 superconductor. The results, summarized below, can be easily generalized for 
any nontrivial pairing state and may be relevant for heavy-fermion superconductors with 
impurities as well. 

Our main results are as follows: (i) A strongly-scattering scalar impurity is a requirement 
for a locahzed, virtual or marginally real, bound state to exist in a li-wave superconductor. 
It is intuitively obvious that any strong enough pair-breaking impurity — magnetic or non- 
magnetic — will induce such a state. Indeed, the low-lying quasiparticle states close to 
the nodes in the energy gap will be strongly influenced even by a non-magnetic impurity 
potential, resulting in a well-defined bound state in the unitary limit. This should be 
compared to the fact that, in s-wave superconductors, both magnetic and resonant non- 



magnetic impurities produce bound states inside the energy gap |10|. (ii) The energy f2' 
and the decay rate Q" of this state are given by 

in 1 



f] = fi' + 2fi" = Ao- ''''^^ 
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where c = cot 6o- We have assumed impurity scattering to be close enough to the unitary 
limit so that the result can be computed to logarithmic accuracy with log(8/7rc) ^1. It 



is only in this limit that the bound state is well-defined. In the unitary limit, defined as 
f^o -^ 7r/2 (c ^ 0), the virtual bound state becomes a marginally bound one at f] ^ with 
Q" /Q' -^ 0. In the opposite case of weak scattering with c ^ 1, the energy of the virtual 
bound state formally approaches Q' ~ Aq and the state is ill-defined because Q" ~ Q' 
(see Fig. 1). The wave function of the bound state is found to decay exponentially, except 
along the directions of the vanishing gap. (iii) For a finite impurity concentration riimp, 
we recover the results obtained earlier for non-conventional superconductors @J^. While 
generally one finds that an impurity band is formed after averaging over impurity positions, 
the Mott criterion for the formation of the impurity band should be modified because of 
highly anisotropic impurity states; we suggest that the metal-insulator transition can be 
observed in the impurity band. 

Single- Impurity Problem. Consider the single scalar impurity problem with H-u^t = 
^kk'a'^'^k ^k'ai whcrc V is the strength of the scalar impurity potential at r = 0, resulting 
in s-wave phase shift ^o- 

The scattering of quasiparticles off the impurity is described by a T- matrix, r(c(j), which 
is independent of wavevector. The Green's function in the presence of an impurity is 
G^. (cu) = Gf{ij)5^j:, + Gf{u)f{u:)Gf{u}), where both Gf{u}) and f (w) are matrices 
in Nambu space. Here [G^ (cj)]"^ = ujtq — A^fi — ^^ts, where ^^ is the quasiparticle en- 
ergy, A^ = Aq cos2(y9 is the gap function with (i^.2_^2 -symmetry, fj {i = 1, 2, 3) are the Pauli 
matrices, and fo is the unit matrix in Nambu particle-hole-spinor space. 

From the previous analysis [@,H,0, it is known that T = TqTq for s-wave scattering and 
a particle-hole symmetric system. Therefore, the only relevant term in the T-matrix takes 
the form To(cu) = Go{uj)/[c^ - Go{ujf], where Go(cu) = 2^ E^Tr G?^(c<;)fo. The virtual 
and bound states in the single-impurity problem are given by the poles of the T-matrix 

c = ±Go(fi), (2) 

which is an implicit equation for fi as a function of c, the strength of impurity scattering. The 
two signs in Eq. (^ are a result of the particle-hole symmetry. Choosing the gap function 



at the Fermi surface so that A{ip) = Aocos2(y9, one finds Go{uj) = (iu/JA{(py — ^^ , 
where the angular brackets denote averaging over the Fermi surface; for simphcity, we take 
{•)fs = /• dip/2n ||rT|. For \uj\ ^ Aq, one finds 



2uj f 4Ao m\ 

In principle, the solution of Eq. (^) is complex, indicating a resonant nature of the quasipar- 
ticle state, better described as a virtual state. This is easily seen form Eq. (|l|), which solves 
Eq. (H) to logarithmic accuracy. However, as c ^ 0, the resonance can be made arbitrarily 
sharp. For c = 0, the virtual state becomes a marginally well-defined state bound to the 
impurity. Exact numerical solution of Eq. (0) as a function of c is shown in Fig. 2. To our 
knowledge, this result has not been claimed before. As c —>■ 1^ , Q' and Q" increase without 
bound so that Q"/Q' —>■ 1~, and the solution becomes unphysical. For c > 1, no solution 
has been found for Q p]2| . 

There are important physical implications of these bound states in a tZ-wave supercon- 
ductor. Consider the most interesting case of unitary impurities in the dilute limit, separated 
by a distance greater than the coherence length ^. Before averaging over impurities, these 
bound states are nearly localized close to the impurity sites (see below) and can substan- 
tially modify the local characteristics of the superconductor: for example, the local density 
of states and the local NMR relaxation rates of atoms close to the impurities. 

Consider a local density of states, defined as Nif^oj) = —^lm.G{r,r;u! + iO^), with 
the total Green's function in the presence of the impurity G{f,f';uj) = G^^''{r — r',uj) + 
G'^'^^f, uj)T(uj)G^^\r', iu), the second term describing the local distortion due to the impurity. 
Using the eigenstates representation of G(r, r'; u) = Z]n('^n(r)V'n(^'))/(^ ~ -^n), "we find two 
terms in the local density of states N{r,uj) = N{uj) + Ni^p{r,u!) ~ J2n{i'ni'^i^n{r))S{^ — En) , 
for well defined states. The first term originates from the bulk quasiparticles, which are 



described by plane-wave eigenstates with E^ = J^i + A|, G^^\f,Lu) = Y.kiu'l/i^ — E^) + 
v'i/{uj+Ej:)], where Uj^ and v^ are the standard Bogoliubov factors. The bulk density of states 
is constant in the system with N{lj)/No = uj/Aq, for u <^ Aq. The second term, A^imp(^, ^) = 



— -lm.[&^\r,uj)T{u!)&^\f,uj)]ii, originates from the virtual quasiparticle state created at 
the impurity: Nimp{r,u) = -^Im^„(V^*jjp,,(r)?/^imp,„(r))/(u;-Eimp,n+«0+). As an important 
example, consider the limit of unitary scattering for which the resonant state is formed at 
Eimp,n = fi — > 0. Because linG^^\f,uj = 0) = —itN(uj = 0) = 0, only the imaginary part of 
the T-matrix contributes to A^imp and the bound-state probability density is found to decay 
as the inverse second power of the distance from the impurity along the nodes of the gap 
function, 

A^imp(r, tu = 0) = Re [G^'^\f, u = 0)]^ oc r'^, (4a) 

and exponentially in the vicinity of the extrema of the gap function, 

iVimp(r,cu = 0) (X [av)/r]-'e-'^/^^^\ (4b) 

where ^(v?) is the angle- dependent coherence length of the superconductor, naturally defined 
as C,{f) = hvF/\'^{f)\- The fact that the impurity state is marginally bound is reflected in 
the logarithmically divergent normalization. This divergence should be cut off at an average 
distance between impurities at any finite concentration. More generally, for an arbitrary 
position of the resonance, taking into account that only one state has been produced with 
^imp,n = n' + in", we find N;^^{f,uj) = lEiF{f- rl)^'/[(cj - n'-f + nf], where we 
have introduced the sum over different impurities, located at r^, and F{f — fl) = {'ipi^pir — 
ri)4'imp{r — fj)) is the probability density of the i-th. impurity state. 

The local variations of the density of states can be probed directly, in principle, by 
scanning-tunneling microscopy. However the NMR experiments on Cu in Zn-doped cuprates 
are quite revealing as well. From Eq. ( |^ ) and below, one concludes immediately that the 
local NMR signal would show two distinct relaxation rates (or even the hierarchy of rates): 
one coming from the Cu sites, far away from the impurities, and another from the sites, 
close to the impurities. The Cu sites near the impurities will be sensitive to the higher local 



density of states and will have a higher relaxation rate at low temperatures [T^. At finite 



impurity concentration (~ 2%), the volume-averaged density of states will have a finite limit 



at cj ^ 0, as follows from Eq. (Pa|). The relaxation rates of Cu atoms close to and away 
from an impurity will, therefore, have the same temperature dependence (TiT)~^ = const, 
but will be of a different magnitude. 

This behaviour has been observed experimentally: Ishida et al. [^ have measured two 
NMR relaxation rates for Cu in Zn-doped YBa2Cu307_5. The second NMR signal with 
higher relaxation was inferred arising from the near-impurity Cu sites. A direct comparison 
of our prediction for local quantities, as probed by NMR, will require a specific model and 
is left for a future publication. 

We would like to contrast our picture of the dilute limit of strongly scattering centers 
to the usual approach of averaging over impurities at finite concentration. If one considers 
averaging over impurities, two NMR relaxation rates, arising from unequivalent sites, cannot 
be resolved; the local inhomogeneous aspect of the localized states is lost after averaging 
over impurity positions. 

For practical purposes the distinction between the impurity bound states and continuum 
in our case is not as well defined as in s-wave superconductors. Any finite temperature 
will produce a finite lifetime for these bound states, and they will be hybridized with the 
continuum of low-energy quasiparticles. 

Finite Concentration of Impurities. Consider the growth of the impurity band with 
finite concentration of strongly scattering impurities. As was mentioned above, scalar (non- 
magnetic) impurities are pair-breakers for any nonconventional superconductor, and they 
substantially change the low-energy spectrum of superconducting quasiparticles. This prob- 
lem has been addressed earlier in great detail (for example, see [0,^). Here we will repeat 
the main steps and give results for the quasiparticle scattering rate and low-energy density 
of states for completeness. 

For finite impurity concentration, the self-consistent Green's function, averaged over 
impurity positions, obeys the Dyson equation {G^{uo))^^ = {G^{ijj))^^ — T,{uj) with 
S(co') = nimpT(ci;) in the single-site approximation. The T- matrix has to be determined 



self-consistently with {Gq^u)) = 2^;^I]^Tr {Gj:{ij))fo. As above [see the remark preceding 
Eq.(^], only Sq = |TrSfo is nonzero. The algebra is straightforward, and for unitary 
scattering yields 

7 ~ y'r2imp(Ao/7riVo), (5) 

where 7 = — IniS(ci; -^ 0) is the scattering rate for low-energy quasiparticles. For cj ^ 7, 
the density of states is determined by impurities and is finite: A^imp(0)/iVo = 27/7rAo. The 
characteristic width of the impurity-dominated region is cj* ^ 7 oc ^n-u^p. 

The origin of the finite density of states is the impurity band, grown from the impurity- 
induced bound states (consider c = 0). Scaling of the impurity bandwidth 7 oc ^yni^ has 
been obtained earlier for the case of paramagnetic impurities in an s-wave superconductor 
0. The fact that 7 oc y^ri^np is obeyed in the case of a d-wave superconductor with scalar 
impurities as well supports our claim that the low-energy states in a disordered d-wave 
superconductor are indeed formed from the bound states at finite concentration. 

Here we shall comment on the condition for the formation of the impurity band with 
extended states, i.e., the Mott criterion. Consider the case of a (3D) s-wave supercon- 
ductor. A magnetic impurity generates a bound state with the impurity wave function 
l^impl ~ exp(— ar), a ~ ^^^ @)0]- For the true conduction band to be formed, the overlap 
between localized states should be large enough. This leads to the Mott criterion for the min- 
imum concentration nj^pO;"^ > 0.2 [|T^. In practice, for conventional superconductors this 
implies a very low critical concentration n? ~ {(^/O^ i^^ being the lattice constant), as the 
coherence length ^ is relatively large ||15| . The situation is qualitatively different in a d-wave 



superconductor. The wave function of impurities has lobes sticking out in the directions 
of vanishing gap, i.e., when cos2ip = [see Eq. (^al)], and the overlap between impurities 
is larger along these directions. The Mott transition in a system of strongly anisotropic 

localized states is an interesting problem and it should be expected that the Mott criterion 

1/2 
should be modified with a somewhat lower constant in the r.h.s. of n^'^ = const and it 



follows that nf oc {a/Ci^ ~ 1%, which is within the current experimental range. A short 
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coherence length of high-Tc superconductors (^ ~ 20A) suggests the possibihty of the ex- 
perimental observation of the metal-insulator transition in the impurity band in a d-wave 
superconductor (assuming that localization effects are small). 

In conclusion, we find that a strongly scattering potential impurity produces a resonant 
or a marginally bound state inside the gap in a rf-wave superconductor. The wave function 
of the impurity bound state is highly anisotropic with 1/r decay along the nodes of the gap 
and exponential with angle-dependent decay range otherwise. These bound states change 
the local density of states Ni^p{r,uj) dramatically, which could be probed experimentally, 
e.g., in NMR. The short coherence length of high-Tc superconductors may lead to a finite 
critical concentration nf within the range of a few percents and allow a direct observation 
of the metal-insulator transition in the impurity band. 
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FIGURES 

FIG.l The relative density of states N{uj)/No for the rf-wave superconductor in the absence 
of impurities (thin sohd hne), the impurity- induced bound state for c — > (thick sohd hne), 
and the virtual bound state for c > (dashed line); A^o is the normal-state density of states 
at the Fermi energy. The finite lifetime {Q" 7^ 0) of the virtual bound state in the d-wave 
superconductor results from the finite density of states A^(u;) oc 00 for small 00 from nodal 
quasiparticles, in contrast to the true bound state in the s-wave superconductor. Exactly 
at u; = 0, the density of quasiparticle states is zero and the virtual bound state becomes 
marginally bound on the edge of the particle-hole continuum. 

FIG. 2 The energy Q = Q'+iQ" of the virtual bound state in the one-impurity problem, given 
by Eq. (0), as a function of impurity strength c: the shown quantities are Q' (solid line), Q" 
(dashed line), and Q" /Q' (chain-dashed line). A spherical Fermi surface and Ag = Aq cos 2^9 
have been assumed. The width Q" of the virtual bound state is always smaller than its 
energy Q' in the neighbourhood of unitary scattering. In contrast, for weak scattering c ~ 1, 
Q" ~ Q' and an impurity-induced virtual state does not exist. The inset shows a comparison 
between the exact result and the asymptotic approximation (dotted lines), as computed to 
logarithmic accuracy by Eq. ([l|) for Q. 
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